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SUMMARY 

Landscape epidemiology and landscape genetics combine advances in molecular techniques, spatial analyses and 
epidemiological models to generate a more real-world understanding of infectious disease dynamics and provide powerful 
new tools for the study of RNA viruses. Using dog rabies as a model we have identified how key questions regarding 
viral spread and persistence can be addressed using a combination of these techniques. In contrast to wildlife rabies, 
investigations into the landscape epidemiology of domestic dog rabies requires more detailed assessment of the role of 
humans in disease spread, including the incorporation of anthropogenic landscape features, human movements and socio- 
cultural factors into spatial models. In particular, identifying and quantifying the influence of anthropogenic features on 
pathogen spread and measuring the permeability of dispersal barriers are important considerations for planning control 
strategies, and may diff'er according to cultural, social and geographical variation across countries or continents. Challenges 
for dog rabies research include the development of metapopulation models and transmission networks using genetic 
information to uncover potential source/sink dynamics and identify the main routes of viral dissemination. Information 
generated from a landscape genetics approach will facilitate spatially strategic control programmes that accommodate for 
heterogeneities in the landscape and therefore utilise resources in the most cost-effective way. This can include the efficient 
placement of vaccine barriers, surveillance points and adaptive management for large-scale control programmes. 

Key words: Canis familiaris, landscape genetics, phylogeography, spatial heterogeneity, transmission dynamics, 
vaccination, zoonosis. 



INTRODUCTION 

Landscape epidemiology is the study of the causes 
and consequences of spatial variation in disease 
incidence or risk across heterogeneous landscapes 
(Ostfeld et al. 2005). Landscape structure affects 
the distribution, abundance and movements of 
host, vector and pathogen populations and therefore 
inherently influences localised interactions between 
infectious and susceptible individuals (Ostfeld et al. 
2005; McCallum, 2008). Revealing the landscape 
factors underlying these interactions calls for an 
interdisciplinary approach that draws on a range of 
techniques across different spatial scales (Manel et al. 
2003; Ostfeld eta/. 2005). Molecular markers provide 
a basis for this by genetically tracking spatial and 
temporal dynamics in pathogen and host populations 
(Biek and Real, 2010). A landscape genetics approach 
to infectious disease therefore encompasses a range 
of analytical tools, including geographic infor- 
mation systems, remote sensing, population genetics, 
phylogenetics and statistical and mathematical 
modelling techniques (Manel et al. 2003). 

* Corresponding author: Kirstyn Brunker, k.brunker.l@ 
research.gla.ac.uk 



RNA viruses represent an ideal group for exploring 
landscape influences on evolutionary trajectories 
due to their characteristically high mutation rates 
and short generation times, which means that 
epidemiological and population genetic processes 
occur on a similar timescale (Drummond et al. 
2003). The accumulation of mutations over time 
and space imprints on the structure of viral genomes 
in a population and can be visualised in data collected 
over months or years, providing a valuable resource 
for the elucidation of ecological and evolutionary 
dynamics. Despite this, the processes that govern 
phylogeographic patterns in viruses are still poorly 
understood (Holmes and Grenfell, 2009), pointing to 
the need for more detailed study into the effect of 
spatial heterogeneity on viral transmission. 

Modern sequence analysis has the power to reveal 
the historical emergence of pathogen variants, distri- 
bution patterns, spillover events, interspecies trans- 
mission and changing selection pressures (Archie 
et al. 2009; Anderson et al. 2010). Advances in 
sequencing technologies have paved the way for a new 
generation of approaches to the study of disease 
dynamics, and next generation sequencing (NGS) 
techniques are being continually refined, improving 
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their resolution, throughput potential, affordability 
and accessibility (Holmes and Grenfell, 2009; 
Metzker, 2010). We envision a paradigm shift to 
whole genome sequencing as the standard technique 
for characterizing RNA viral evolution on small 
spatio-temporal scales, providing greater discrimi- 
nation between genotypes and finer resolution in 
population structure (Holmes and Grenfell, 2009). 
Given the progress in sequencing technology and 
parallel advances in spatial analytical tools, it is an 
exciting time to study the landscape epidemiology of 
RNA viruses from a population genetic perspective. 

In this review, we focus on rabies virus (RABV), 
which presents an excellent model system to illustrate 
the challenges and prospects of such an approach. 
RABV is a single stranded, negative-sense RNA virus 
belonging to the Lyssavirus genus (Family: 
Rhabdoviridae) (Dietzschold et al. 2005). It is 
globally distributed and has the ability to infect all 
mammals, but typically exists in endemic foci as a 
reservoir host-specific variant with occasional spil- 
lover to other species (Rupprecht et al. 2002). Domestic 
dogs, Canis familiaris, are the principal reservoir of 
RABV, responsible for 99% of the estimated 55,000 
human deaths due to rabies that occur mainly in Asia 
and Africa every year (Knobel et al. 2005). 

Rabies has proven a remarkably valuable system 
for exploring the effect of landscape processes on 
host/pathogen interactions. However, most studies 
have focused on wildlife rabies due to the higher 
quality of surveillance data and availability of 
resources in areas with major endemic wildlife 
foci, e.g. raccoon rabies in eastern North America 
or fox rabies in Western Europe (Bourhy et al. 1999; 
Smith et al. 2002; Holmes 2004; Biek et al. 2007; 
Wheeler and Waller 2008; Szanto et al. 2011). The 
wealth of research into wildlife rabies provides a basis 
for comparison with the domestic dog foci that exist 
in Africa, Asia and parts of Latin America, which, 
despite their much greater public health burden, have 
been less well studied. 

Domestic dogs are inherently tied to human 
populations, and various aspects of human ecology, 
including distribution, habitation and movement 
patterns, or cultural practices, will directly influence 
rabies spread in dog populations. Settlements can be 
considered dog 'habitat', and dog densities have been 
predicted on the basis of human demographics and 
human geography (Butler, 2000; Knobel et al. 2008). 
A seemingly ubiquitous feature of countries with 
persistent dog rabies foci is the free-roaming nature 
of these populations, often referred to as 'neighbour- 
hood' dogs. Dogs in Africa, Asia and parts of Latin 
America, where canine rabies is endemic, are rarely 
restricted by leashing or enclosures and their role 
as domestic animals varies e.g. watch dogs, trade, 
companion animals. The free movement of dogs 
would thus be expected to contribute to local rabies 
transmission, potentially resembling the known 



features of wildlife rabies and illustrating the complex 
interplay between anthropogenic and natural drivers 
of disease spread in this system. A key challenge is to 
uncover the extent to which natural constraints to 
rabies flux hold in a host-pathogen system with 
greater human-mediated dispersal; and what affect 
this has on phylogeographic signatures. 

This review aims to synthesise our current 
understanding of rabies landscape epidemiology 
derived from genetic data across different spatio- 
temporal scales (see Table 1 for a list of relevant 
studies and associated analytical methods). Using 
the extensive work on wildlife rabies as a backdrop, 
we seek to identify commonalities as well as 
fundamental differences characterizing the dyn- 
amics of RABV in domestic dog populations. We 
will argue that the specific 'landscape' supporting 
the sustained transmission of the virus in dogs is 
determined by a complex mixture of physical and 
human geography and that a quantitative under- 
standing of these landscapes will be essential for 
rabies control and eradication. In addition, control 
activities are themselves predicted to change mol- 
ecular epidemiological trajectories, creating interest- 
ing opportunities for adaptive management in 
rabies. Throughout the review we highlight how 
novel tools and technologies are being used to tackle 
these problems and identify key areas in which such 
approaches have future potential. 

RABIES LIFE CYCLE AND PHYLOGEOGRAPHY 

Transmission of rabies virus occurs predominantly 
through the bite of an infected animal, which 
inoculates virus-laden saliva into the subcutaneous 
and muscle tissue of a susceptible host (Dietzschold 
et al. 2005). Once inoculated, the virus enters neurons 
and migrates to the central nervous system, before 
spreading to other organs, including the salivary 
glands where large amounts of infectious virions are 
shed into the saliva for further transmission 
(Dietzschold et al. 2005). The incubation period is 
highly variable with a mean of 22 days in naturally 
infected dogs (Hampson et al. 2009), but may extend 
to months or years due in part to localised replication 
(Hanlon et al. 2007). In contrast, the infectious 
period is very short, around 3 days, and very rarely 
exceeds 10 days. The short infectious period severely 
restricts the spatial scale over which an infectious 
individual can transmit the virus, while the longer 
and more variable incubation time may permit the 
active or passive movement of infected individuals 
over larger distances, which, as detailed below, helps 
to explain phylogeographic patterns. 

On a global scale, dog RABV forms six major 
genetic clades, five of which are associated with 
particular geographic regions: Africa 2, Africa 3, 
Arctic-related, Asian, and the Indian subcontinent 
(Bourhy et al. 2008). Grouping of these clades largely 



Table 1 . Studies that have examined sources of spatial heterogeneity in dog rabies dynamics at a landscape scale; Gen = genetic data, Epi = epidemiological data. 

ML = Maximum Likelihood 



Scale 


Sources of 
spatial 

heterogeneity 


Data 


Analytical methods 


Key points/Summary 


Reference 


Local ( ~ 5 yr 


Individual host 


Epi 


Reconstruction of epidemic trees, outbreak 


Contact tracing data used to generate 


Hampson 


period) 


heterogeneity 




simulations, construction of transmission 


robust estimates of epidemiological 


et al. 2009 








networks based on a spatial infection 


parameters 










kernel and generation intervals estimated from 












epidemiological data using ML 






Local ( ~ 5 yr 


Spatial configuration of 


Epi 


Patch-occupancy models 


Uses dog bite incidence records to explore 


Beyer et al. 


period) 


populations 






metapopulation dynamics 


2011 


Regional 


Socio-economic drivers 


Gen 


Bayesian inference of phylogeny, molecular 


Attributes phylogeographic patterns to economic 


Carnieli et al. 


(20+yrs) 






clocks and demography 


growth and migration patterns of humans 


2011 


Regional 


Translocation, road 


Gen-Epi 


Phylogenetic inference 


Patterns suggest importance of wildlife hosts in this 


Coetzee and 


(10+ yrs) 


networks, wildlife hosts 






system. Spread of rabies coincided with major 


Nel, 2007 










highways, and indicated translocation events. 




Regional 


Translocation 


Gen 


Antigenic characterization through monoclonal 


Identified the regional source of a rabid dog case 


David et al. 


(2 yrs) 






antibody profiling and molecular sequence 


through forensic epidemiological tracing 


2004 








comparison 






Regional 


Reservoir hosts, cultural 


Gen 


Bayesian phylogenetics; parsimonious 


Used statistical parsimony to construct 


Lembo et al. 


(10 yrs) 


drivers 




construction of transmission networks 


most likely transmission networks. Also 


2007 










discovered the potential influence of social 












drivers on rabies phylogeographic structure 












(pastoralist vs. non-pastoralist community 












structure). 




Local and 


Reservoir hosts. 


Gen-Epi 


Bayesian phylogenetics ; generation of 


Most likely reservoir hosts (domestic dogs) 


Lembo et al. 


Regional 


transmission clusters in 




epidemic trees based on probabilities of 


inferred. Found evidence for short-lived 


2008 


(10+ yrs) 


wildlife 




links between possible progenitors and 


chains of transmission in wildlife. 










suspected cases weighted by spatio-temporal 


compared to self-sustaining transmission 










proximity, ML estimation of spatial infection 


in domestic dogs 










kernel and generation interval distribution 






Regional 


Road networks. 


Epi 


Spatial analysis of reported rabies using GIS. 


Distribution of cases followed road network 


Tenzin et al. 


(1 yr) 


population density 




Directional spread of rabies based on mean 


and towns with high human density and high 


2010 








centre of cases and a standard deviational ellipse 


numbers of free-roaming dogs 










weighted by date of cases 






Country 


Translocation, cross 


Gen-Epi 


Bayesian inference of phylogeny, molecular clock, 


Frequent incursions across country boundary; 


Hayman et al. 


(3 yrs) 


border incursions 




demography and phylogeographic diffusion 


first molecular evidence for a long distance 


2011 








based on discrete spatial states 


translocation of a rabies sub-lineage in Africa 
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reflects major barriers such as oceans, large mountain 
ranges and deserts, or historical colonisation events. 
The Arctic-related viruses are relatively well dis- 
tributed, reflecting the lack of barriers in the far 
north, and a sixth clade has a cosmopolitan distri- 
bution, reflecting historical waves of human 
migrations and colonisations (Smith et al. 1992; 
Bourhy et al. 2008). Such efficient dissemination 
is likely to be a result of the virus's relatively long and 
variable incubation period occasionally facilitating 
long-distance transport of infected dogs. 

The strong phylogeographic structure exhibited 
by most dog rabies clades is mirrored in wildlife 
rabies, where the same pattern is evident even at 
relatively small spatial scales (Biek et al. 2007) and 
may be explained by a 'surfing mutation' model 
(Excoffier and Ray, 2008). According to this model, 
new lineages arising during initial colonisation are 
able to reach high frequencies driven by an epidemic 
wave, whereas subsequent lineages do not benefit 
from such conditions and fail to infiltrate the 
dominating lineage clusters. Initial invasion events 
can therefore markedly influence the phylogeo- 
graphic structure of rabies and, at least in wildlife, 
this structure can remain intact for decades (Biek 
et al. 2007; Szanto et al. 2011). However, with 
ongoing movement of infected individuals and the 
immigration of new lineages into an area, this 
structure would be expected to erode over time. 
How environmental variation affects these processes 
of emergence and subsequent erosion of phylogeo- 
graphic structure are questions central to defining the 
landscape epidemiology of RABV and we will return 
to them throughout this review. 

LANDSCAPE LEVEL EFFECTS ON RABIES 
DYNAMICS 

Since the transmission of rabies is entirely dependent 
on host movements initiating contact between in- 
fectious and susceptible individuals, the landscape 
that they occupy and disperse within heavily influ- 
ences the ability of the virus to infect new hosts. 
Indeed, wildlife rabies epidemics tend to spread as 
irregular waves that differ in velocity according to 
heterogeneities in the landscape (Russell et al. 2006). 
These 'heterogeneities' include natural features such 
as mountain ranges, water bodies and deserts, but also 
anthropogenic features including roads and vaccine 
corridors. The influence of spatial heterogeneity on 
rabies spread can be broken down into three aspects, 
each of which are discussed in detail in the following 
sections: (1) host movements: natural versus human- 
mediated; (2) landscape attributes influencing rabies 
spread; and (3) population level effects. 

The role of human vs. natural dispersal 

Throughout history, humans moving animals has 
repeatedly led to the emergence and spread of rabies 
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Association index (AI): test statistic used as a means of quantifying phylogeny-trait associations i.e. given a discrete 
character for each tip of a phylogenetic tree, are more closely related taxa likely to share the same trait values than is 
expected by chance alone. Avoids the issue of lack of independence due to shared ancestry (Wang et al. 2001). 

Metapopulation: concept describing the persistence of a species in a spatially heterogeneous environment as a balance 
between colonisation and extinction in loosely coupled subpopulations or 'patches' with different within and between- 
patch dynamics, which can be applied to infectious diseases (Grenfell and Harwood, 1997) 

Most recent common ancestor (MRCA): most recent individual from which all taxa in the group are directly 
descended. 

Oral rabies vaccination (ORV): distribution of oral rabies vaccine baits as a strategy to control the spread of wildlife 
rabies. Used as a control measure for wildlife rabies in North America, and Europe. 

Phylogeography: combines phylogenetics and biogeography to describe the contemporary pattern of an organism's 
geographic spread according to gene genealogies. 

Ro (basic reproductive number): the average number of secondary cases derived from a single infectious individual 
in an entirely susceptible population (Anderson and May, 1991). 

Superspreader: an individual causing a disproportionately high number of secondary cases, compared to the mean 
(represented by Rq), depicted in the long tail of a frequency distribution of secondary cases (Lloyd-Smith et al. 2005). 

Surfing mutation model: genetic variants at the front of an advancing wave of infection are swept to high frequencies 
during an epidemic peak, resulting in long-term dominance of colonising genetic lineages (Excoffier and Ray, 2008, 
Klopfstein et al. 2006). 



Box 2. Key research questions for dog rabies at the landscape scale 



1. To what degree do the demography and ecology of dogs, as compared to humans, determine the 
dynamics of rabies transmission? 

As dogs are inherently tied to humans, we might expect the patterns dictated by local dog movements to be 
confounded by human-mediated long-distance movements or anthropogenic features that facilitate connection of 
sub-populations. 

2. Which aspects of human geography best predict landscape permeability to dog rabies? Are these 
predictors consistent across different areas and continents? 

Although evidence suggests the influence of human geography on dog rabies dispersal, we have yet to uncover the 
best predictors of this form of spread. Specifically, a quantifiable method of characterizing these landscape features 
would assist with the creation of guidelines for using dispersal barriers to aid control. 

3. How does endemic rabies compare to epidemic spread? 

Most studies have focused on outbreak situations in wildlife, but it is not clear how rabies dynamics change in 
systems where the pathogen has been circulating for centuries, as is the case for dog rabies. How long is 
phylogeographic structure maintained over time and do initial invasion pathways predict connectivity in the 
endemic state? 

4. Does dog rabies persist in metapopulations? 

Dog rabies may best be described as a series of sub-populations with varying inter- and intra-patch dynamics. 
Exploring the concept of a metapopulation dynamic for dog rabies presents a relatively unexplored area for future 
research and may uncover the mechanisms that allow pathogen persistence at a local scale. 

5. How can information about landscape heterogeneity and genetic structure be incorporated into more 
efficient control programmes? 

Information gathered from spatially explicit landscape models will allow targeted control measures using resources 
cost-effectively. Moreover, the implementation of vaccine barriers and other forms of control change the complexity 
of the landscape, potentially altering disease dynamics. Powerful genetic techniques may elucidate the effect of these 
new landscape heterogeneities, facilitating adaptive management. 



in susceptible populations, either through deliberate 
translocation (e.g. the relocation of animals, some of 
which may be incubating virus (Nettles et al. 1979)) 
or inadvertent movement (e.g. raccoons on garbage 
trucks (Wilson et al. 1997)). 

Although human-mediated long-distance move- 
ment was found to be a significant feature of raccoon 
rabies spread (Smith et al. 2002), it was still rare 
compared to natural dispersal, and stochastic in its 
occurrence. In contrast, human-mediated dispersal 
of dog rabies appears more widespread and 



potentially predictable due to possible links with 
human activities and migration patterns. The 
European colonisation of Africa is thought to have 
led to the introduction and subsequent expansion of 
canine rabies across the continent (Bourhy et al. 2008; 
Talbi et al. 2009; Lemey et al. 2010); increasingly 
frequent records of re-introductions into countries 
where rabies has been eliminated have been recorded 
(Zanoni and Breitenmoser, 2003; Weiss et al. 2009; 
Gautret et al. 2011;), and canine rabies has recently 
emerged on several previously rabies-free islands in 
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Indonesia via fishermen importing incubating dogs 
(Windiyaningsih et al. 2004; Susilawathi et al. 2012). 
Even on smaller scales, phylogeographic patterns 
frequently imply translocation events that increase 
the level of mixing between lineages and may obscure 
patterns of lineage clustering (David et al. 2004; 
Talbi et al. 2010; Hayman et al. 2011). 

Although often inferred, it is difficult to quantify 
the influence of human-mediated movements on the 
diffusion of RABV. Though translocations may be 
frequent, the likelihood of establishment and spread 
is much lower (Smith et al. 2005). But, unlike for wild 
animals, where translocation is risky and often 
unsuccessful, human-mediated displacement from a 
dog's original home range confers fewer risks. Dog 
translocations are therefore likely to be more success- 
ful in initiating new disease foci, and with potentially 
higher consequences, introducing disease into pre- 
viously disease-free areas. Given that long-distance 
movements instigated by humans are often in 
response to social drivers, the genetic structure of 
dog rabies is expected to reflect variation in cultures 
and socioeconomic factors. For example, rural 
workers in Thailand often relocate to urban areas 
where work is more readily available in the off- 
growing season, taking their dogs, and sometimes 
rabies, with them (Denduangboripant et al. 2005). 
Similarly, the resurgence of canine rabies in 
KwaZulu-Natal in the 1970s was linked to refugee 
movements from Mozambique (Swanepoel et al. 
1993; Cleaveland, 1998) and phylogenetic patterns 
of RABV in parts of northern Tanzania can be 
tentatively explained by the movement of nomadic 
Maasai pastoralists with their dogs (Lembo et al. 
2007). Counterproductively, dog-owners have been 
reported to move their dogs to avoid the threat 
of culling during attempts to control rabies, which 
may explain more rapid spread of rabies than would 
be predicted by dog movement alone. Importantly, 
many of these aspects of human geography are 
quantifiable, making it possible in principle to 
generate testable hypotheses about their role in 
determining rabies phylogeographic structure. 

Talbi et al. (2010) used spatial simulations to show 
that the observed patterns of spread in North Africa 
could only be explained by the inclusion of long- 
distance translocation events as opposed to natural 
dog movements alone. The wave-like patterns that 
are a signature of wildlife rabies and determined by 
natural host movements are less evident in domestic 
dog rabies, and it is unknown to what extent natural 
dog movements vs. human-mediated movements 
determine the observed phylogeographic patterns. 
A recent study in a Kenyan rangeland found 
that healthy domestic dogs rarely moved more than 
50 metres from their home bases, and the maximum 
distance recorded was 3-2 km (Woodroffe and 
Donnelly, 2011). In contrast, dispersal distances in 
wildlife hosts tend to be at least an order of magnitude 



higher (Cullingham et al. 2008). However, rabies 
infection leads to behavioural changes, which may 
alter movement patterns; in northern Tanzania, most 
contacts with rabid dogs occurred within a kilometre 
of an animal's homestead (mean ~0-88km), but a 
small proportion of rabid dogs ran over 1 5 km while 
infectious (Hampson et al. 2009). Dog movements 
also vary depending on the societal context and 
livelihoods of the communities to which they belong 
(Woodroffe and Donnelly, 2011). For their dispersal 
simulations, Talbi et al. (2010) used measures of 
natural dog movements based on a specific locality in 
Tanzania (Hampson et al. 2009) but it remains 
unclear whether such data are transferrable among 
different geographic localities. Research on domestic 
dog ecology and movement is lacking and detailed, 
location-specific data are required in order to reliably 
test hypotheses for alternative drivers of viral diffu- 
sion. 

Landscape attributes influencing rabies spread 

As a parasite transmitted through direct host-to-host 
contact, spatially-defined genetic discontinuities 
in RABV populations may indicate a barrier to 
host contact (Biek and Real, 2010). We consider a 
'dispersal barrier' to be any spatial feature of the 
landscape that impedes the gene flow of a pathogen. 
Spatially explicit models incorporating landscape 
features can help uncover barriers to gene flow 
based on slower progression than expected over a 
homogeneous area. However, the interpretation 
of genetic discontinuities requires caution as phylo- 
geographic patterns may arise as a result of historical 
colonisation events and irrespective of physical 
barriers (Real and Biek, 2007; Talbi et al. 2009). 

Smith et al. (2002) demonstrated a seven-fold 
reduction in rates of raccoon rabies spread in North 
America due to the presence of rivers and forest 
cover. Similarly, the Vistula River separated distinct 
clusters of a red fox variant and raccoon dog/fox 
variant in Europe (Bourhy et al. 1999). Clearly, large 
water bodies impede host movement, but have 
differential success as barriers according to additional 
factors such as habitat suitability (Smith et al. 2005; 
Cullingham et al. 2009), physical geography (Rees 
et al. 2009), and width and flow rates (Bourhy et al. 
1999). In an attempt to assess the differential 
permeability of barriers, Rees et al. (2008) used 
genetic simulation modelling to assess the effect of 
the Niagara River on rabies spread from New York 
State to Ontario. Comparing genetic population 
structure derived from field data, with simulated 
population expansion scenarios they ascertained 
a 50% barrier effect. In contrast, homogeneous 
landscapes that lack environmental barriers are 
particularly vulnerable to the rapid expansion of an 
introduced pathogen, as demonstrated by the in- 
creased speed with which raccoon rabies is predicted 



Landscape epidemiology and genetics of rabies in domestic dogs 



1905 



B 




Hull 



11 .4 dogs/km |l 




Fig. 1 . Varying spatial complexity in areas with endemic dog rabies as a result of increasing dog population density, 
A) low density: Ngorongoro District, Tanzania; B) medium density: Serengeti District, Tanzania, C) high density: 
Hermosillo, Mexico. Red circles highlight settlements in rural areas. Maps obtained using Google Earth (littp://earth. 
google.com). 



to cross central Ohio, which lacks major natural 
barriers, compared to neighbouring states (Russell 
et al. 2005). With the exception of major landscape 
features inferred as barriers to dog rabies on a global 



scale (Bourhy et al. 2008), we know relatively little 
about barrier effects at smaller scales. Determining 
whether natural barriers prevent the dissemination of 
dog RABV or if ties to human ecology negate their 
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effect is an important research question with signifi- 
cant implications for control. 

Pohtical borders are another potential form of 
barrier that may play a role in the containment of dog 
rabies, but have no evident impact on wildlife rabies 
dispersal. For example, distinct and almost mono- 
phyletic RABV groups associated with North African 
countries indicate restricted movement across geo- 
political boundaries (Talbi et al. 2010). This well- 
defined population structure at a regional scale 
contrasts with a relatively fluid dissemination within 
countries. On the one hand, this suggests that 
country level vaccination programmes should have 
a good chance of eliminating dog rabies even in 
contiguous landscapes. However, these findings are 
not generalisable, with epidemiological analyses 
of RABV in Ghana highlighting frequent cross- 
border incursions (Hayman et al. 2011). In addition, 
time-series analyses indicate large-scale synchronous 
dynamics of rabies across multiple countries in 
eastern and southern Africa, possibly due to a 
combination of human or wildlife-mediated long- 
distance dispersal and a lack of sustained control 
programmes (Hampson et al. 2007). Understanding 
the circumstances whereby political boundaries act 
as dispersal barriers should provide guidance for 
whether control programmes require regional co- 
operation or can be sustained at a national level with 
appropriate border controls. As a first step, it would 
be useful to compare phylogeographic structure of 
dog RABV among different parts of the world and 
across hierarchical spatial scales. Objective measures 
for such a comparison can be obtained by quantifying 
the degree to which sequences cluster on a phylogeny 
according to their geographic location. Several 
statistics are available for this (Parker et al. 2008), of 
which the association index (Al) (Wang et al. 2001) 
has been found to be of particular utility. This 
measure of phylogeny-trait correlation may 
provide an initial descriptive analysis of consistency 
or variability in the hierarchical genetic structure 
of dog rabies between continents or countries that 
could determine the generalisation of results across 
different systems. For example, we might hypoth- 
esise that regions with large population densities and 
a high level of transport infrastructure, e.g. parts of 
Asia, show less distinctive phylogeographic patterns 
(a low Al value) compared to more sparsely 
populated and less developed landscapes like those 
in sub-Saharan Africa (high Al value). 

Quantifying barrier effects, to predict the like- 
lihood of incursions or to exploit them for control 
programmes, is an important area yet to be tackled for 
dog rabies. The impact of anthropogenic landscape 
features on spread, at least in a qualitative sense, has 
been noted in several countries. For example, Tenzin 
et al. (2010) mapped the spread of rabies in Bhutan 
showing a strong visual pattern alongside road 
networks and towns with high dog-to-human ratios; 



and phylogeographic patterns in north Mexico and 
Thailand match the distribution of major migration 
routes (de Mattos et al. 1999; Denduangboripant 
et al. 2005). It can be hypothesised that features 
important to wildlife rabies become less significant 
to dog rabies as anthropogenic features take over the 
landscape and mediate the effect of barriers, i.e. 
bridges, roads and transportation make it possible for 
dogs to circumvent natural barriers and support a 
relatively fluid dissemination of rabies across natu- 
rally heterogeneous landscapes. These hypotheses 
remain largely untested and thus present a fruitful 
area for landscape genetics investigations. 

Uncovering those features that facilitate viral 
spread is an equally important aspect of assessing 
the landscape. As previously discussed, anthropo- 
genic effects may facilitate the transmission of rabies 
across a larger scale than natural movements 
alone allow. The impact of human-mediated dog 
transport may by explored by phylogenetic analysis, 
with translocations implicated by the presence of 
a cluster-specific variant in a distant locality (see 
David et al. (2004) and Cohen et al. (2007), Table 1), 
essentially an assignment approach (Paetkau et al. 
1995). While this method is not a definitive measure 
of a translocation, it can identify the most likely 
dispersal scenario. Talbi et al. (2010) used Bayesian 
phylogeographic diffusion models (Lemey et al. 
2009) in an attempt to quantify the importance of 
various anthropogenic predictors on the observed 
spread of rabies cases in North Africa (see Table 1). 
Model fitting indicated that dispersal patterns among 
towns were best explained by road distance, consist- 
ent with the anticipated role of human movement. 
Interestingly, further refinements to calculating 
distance matrices, such as an accessibility index 
based on road type and travel time, received only 
limited model support and road distances were 
only a marginally better predictor than great-circle 
distances. This may indicate that RABV dispersal 
follows a rather homogeneous spatial diffusion 
process, without any particular effect of landscape 
heterogeneity. More likely however, it means that 
more useful geographic predictors of human- 
mediated dispersal have yet to be found. While this 
was not possible in the Talbi et al. (2010) study, their 
method can accommodate a wide range of geographi- 
cal and environmental predictors and thus provides a 
promising general framework for examining land- 
scape effects in future data sets. 



Population level effects and metapopulation dynamics 

Infectious disease dynamics are often described in 
terms of a metapopulation, with host populations 
divided into smaller, spatially structured sub- 
populations that exist with different inter- and 
intra-patch dynamics (Grenfell and Harwood, 



Landscape epidemiology and genetics of rabies in domestic dogs 



1907 



1997). Heterogeneities in the spatial configuration of 
host populations are critical to understanding the 
persistence of endemic pathogens (Hagenaars et al. 
2004) and incorporating this social/spatial structure is 
particularly important for developing mechanistic 
models of acute-acting infections like rabies (Cross 
et al. 2005). In the case of dog rabies, human 
settlements (village/town/city) can be considered 
habitat patches that vary in host density and 
connectivity (Fig. 1). For example, rural areas 
typically exist as an array of villages connected to 
larger towns and cities by major roads. 
Fragmentation and low connectivity, i.e. long dis- 
tances between settlements, natural barriers and lack 
of transportation networks, likely restrict the ability 
of a dog to move across a landscape and hence limit 
rabies spread. Hypothetically, limited contact be- 
tween dog populations and correspondingly strong 
spatial structure of rabies would be expected in the 
least developed regions where travel is most limited. 

Local within-patch dynamics are potentially a key 
component contributing to the persistence and 
maintenance of disease at larger scales. Therefore, it 
is important to quantify the relative contribution of 
different patches to overall disease persistence in a 
system. This might involve characterizing patches 
based on dog density, turnover and growth rates, 
levels of ownership and vaccination, alongside 
measures of connectivity. For example, dog turnover 
rates can be extremely high e.g. the Machakos region 
in Kenya has a dog population estimated to grow by 
9% per annum, making rabies control difficult in this 
area (Kitala et al. 2001). Areas such as these may 
possibly be identified as hotspots or source patches 
where rabies has a high chance of being maintained 
and spread from. 

A potentially useful tool for quantifying how 
rabies dynamics differ among local patches is the 
basic reproductive number, Rq. Theory would suggest 
that for directly transmitted diseases such as rabies, 
Rq should increase with host density. However, 
empirical evidence for density-dependent dynamics 
is equivocal, with no detectable differences in Rq 
among dog populations with varying densities 
around the world (Hampson et al. 2009). While it is 
unclear why Rq is so insensitive to differences in dog 
density, e.g. Mexico ( ~ 100 dogs/km^ (Eng et al. 
1993)) versus rural Africa (~ 10 dogs/km^ (Lembo 
et al. 2008)), this finding is more consistent with a 
frequency-dependent mode of rabies transmission. 
For both density- and frequency-dependent disease 
dynamics stochastic extinctions are expected once 
density drops below a certain level, meaning that low- 
density patches may act as barriers to spread. Indeed, 
rabies has been observed to spread but appears 
less able to persist in low density dog populations, 
e.g. Ngorongoro in Tanzania (average of 4-2 dogs/ 
km^ (Lembo et al. 2008)), suggesting that neighbou- 
ring source populations are required for rabies 



maintenance. Cross et al. (2007) discuss the incor- 
poration of heterogeneities in population structure by 
expanding on the utility of Rq as a measure of disease 
transmission. Rq assumes that the population is 
evenly mixed, but the hierarchical nature of disease 
invasion requires stochastic models that incorporate 
both within-patch transmission (Rq) and the factors 
contributing to persistence, i.e. between patch 
transmission - namely the recruitment of suscepti- 
bles, group size and the infectious period. These 
questions of potential metapopulation structure, 
patch variability and source-sink dynamics 
(Pulliam, 1988) represent a particularly pertinent 
and rich area for future studies of RABV which has 
yet to be explored using genetic approaches and are 
vital to inform efficient control strategies. 

Host movement rates and connectivity between 
sub-populations are crucial predictors of disease 
invasion (Cross et al. 2005; Beyer et al. 2011). 
Epidemiological data available at a localised scale, 
e.g. the incidence of dog bites, can provide a simple 
spatio-temporal measure for metapopulation analy- 
sis, as shown by Beyer et al. (201 1). Patch-occupancy 
models identified two metrics of connectivity, the 
distance between neighbouring villages and the size 
of villages receiving infection, as significant factors 
facilitating the transmission of disease. From a 
population genetics perspective, coalescent analysis 
of sequence data could provide the means to identify 
metapopulation dynamics that may not be discern- 
able using epidemiological data alone. Such tech- 
niques have been successfully applied, for example, 
to understand the persistence of influenza A virus, 
with molecular clock-based estimates of divergence 
times of most recent common ancestors (MRCA) and 
demonstrating persistence due to dynamic migration 
patterns rather than source-sink dynamics (Russell 
et al. 2008; Bedford et al. 2010). 

More generally, the fine-grained spatial genetic 
structure evident in wildlife rabies that is lacking for 
dog rabies may be due to differences in time frame. 
The best-studied wildlife rabies dynamics come from 
invasions that originated a few decades ago (Bourhy 
et al. 1999; Real et al. 2005a; Biek et al. 2007), 
compared to endemic foci in domestic dogs, thought 
to have persisted for centuries (Bourhy et al. 2008). 
Factors influencing persistent endemic cycles are 
likely to differ from those determining epidemic 
expansions, with viral population turnover, increased 
mixture of lineages and the greater role of human- 
mediated movement resulting in less clearly defined 
fine-grain viral structure. Again, metapopulation 
models that capture population connectivity at 
appropriate spatial and temporal scales may shed 
light on the genetic structure of endemic foci. 
Characterizing the landscape connectivity between 
rabies endemic areas (or individual rabies cases, given 
the resolution of spatial data) may be facilitated using 
techniques applied in conservation genetics that have 
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yet to be exploited for disease dynamics. For 
example, the program Fractionnator, (http://www. 
unil.ch/biomapper/frictionnator/frictionnator.html), 
provides a 'strip statistic' quantifying the effect 
of landscape features on the relatedness of all 
possible pairs of individuals sampled. This entails 
defining grid cells within a strip across the sampled 
landscape and assessing the abundance of landscape 
features within each strip. In addition, advances 
in the development of spatial analysis software 
such as python-based customised GIS tools (e.g. 
Etherington, 2001) that allow the visualisation and 
measurement of genetic relatedness and landscape 
connectivity based on least cost path analysis, R 
packages such as Maptools (Lewin-Koh et al. 2012) 
and alternative open-source software such as 
PASSaGE 2 (Rosenberg and Anderson, 2011) 
provide utilities for measuring connectivity between 
patches, quantifying patterns in spatial data and 
potentially identifying barriers to dispersal. 

Elucidating how heterogeneity at the host level 
contributes to dynamics at the patch level is an 
important aspect of understanding the mechanisms 
underlying rabies persistence and spread. Differential 
within-patch dynamics resulting from individual 
heterogeneity including variation in biting fre- 
quency, host population size, the structure of a 
settlement, may confer properties that promote the 
effectiveness of a patch as a source of infection. 
Tracing transmission pathways may provide a further 
means of elucidating source populations that initiate 
chains of infection, particularly at very fine-scale 
resolution. Lembo et al. (2007) (Table 1) used 
parsimonious transmission networks to infer reser- 
voir host dynamics and patterns of interspecific 
transmission, but this technique could also be used 
more generally to uncover dispersal patterns, e.g. if a 
particular patch repeatedly acts as a source of 
infection. The construction of transmission trees 
from epidemiological data on the timing and locality 
of rabies cases is another useful approach for 
identifying transmission pathways and can be sup- 
plemented by contact tracing data (Lembo et al. 
2008). Contact tracing of rabies is aided by the 
memorable nature of rabies bites and therefore 
presents a unique situation where the collection of 
contact data is achievable. However, building trans- 
mission trees in this way becomes increasingly 
unreliable as larger numbers of missing links are 
inferred, and these models may also simplistically 
predict unlikely transmission events if knowledge of 
the landscape is not incorporated into algorithms e.g. 
predictions of transmission across major landscape 
barriers. Genetic data could be used to quantify 
mutation rates and build independent time-scaled 
transmission networks against which epidemiologi- 
cally constructed transmission trees could be cali- 
brated. Ultimately, transmission trees will be most 
effectively constructed with the combined use of 



genetic and epidemiological data, and appropriate 
datasets are increasingly becoming available for 
rabies, but further computational and statistical 
advances are required. 

Host transmission is typically heterogeneous, often 
described by the 80/20 paradigm in which individ- 
uals differentially carry and transmit pathogens i.e. 
80% of infections are carried by 20% of the population 
(Woolhouse et al. 1997). In the most extreme case this 
heterogeneity can exist in the form of superspreaders 
(Lloyd-Smith et al. 2005). There have been no in- 
depth studies of superspreading in domestic dogs 
carrying rabies and the contribution of supersprea- 
ders to rabies spread and persistence is controversial 
(e.g. Hampson et al. 2007 vs. Talbi et al. 2009), but 
large variations in biting frequency have been 
observed (Hampson et al. 2009). An important 
distinction must be made between superspreading 
in a restricted area, e.g. the same village, and a 'spatial 
superspreader' that has dispersed infection over large 
distances and multiple settlements (see Fig. 2 for 
visual example). Although the number of secondary 
cases caused by such individuals may be equivalent, 
spatial superspreaders are more important from an 
epidemiological perspective if they are responsible 
for connecting sub-populations. Tracing resulting 
infections to their source through genetic sequence 
data may identify hotspots where certain patch 
features promote greater dispersal ability i.e. spatial 
superspreading events (as suggested by Cross et al. 
2007). Ideally, this would utilise retrospective genetic 
tracing of networks of infection, similar to that used 
for identifying source infections in foot-and-mouth 
disease (FMD) (Cottam et al. 2008). Statistically 
robust methods are yet to be conceived, but account- 
ing for missing links in transmission chains is an 
important area for future development. As a major 
advantage compared to FMD, contact that may have 
lead to transmission is much easier to define for 
rabies. Contact tracing data, if available, could 
therefore complement genetic data to offer a more 
comprehensive view of how individual heterogeneity 
can impact on dynamics. 

INTEGRATING LANDSCAPE EPIDEMIOLOGY INTO 
RABIES CONTROL 

Landscape genetics has undoubtedly generated fun- 
damental insights into the dynamics of rabies but has 
arguably yet to make major contributions to its 
control. The preceding discussion illustrated poten- 
tial avenues for exploration and here we aim to more 
clearly define a landscape genetics research agenda 
that could directly benefit the planning or implemen- 
tation of programmes that aim to control or eliminate 
rabies. 

Considerable progress has been made in the 
development of blueprints and operational toolkits 
for rabies control (see "Blueprint for rabies 
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Fig. 2. Dispersal of bites from superspreading dogs resulting in rabies transmission in an area of the Serengeti District 
in Tanzania. Roads and rivers are shown to highlight the potential influence of landscape features on the dispersal of 
rabies- tentative observations indicate that superspreader progeny appear to cluster alongside roads and movement may 
be restricted by the presence of rivers (but other landscape features not shown may also be responsible for influencing 
dispersal patterns). Two potential types of superspreader are also highlighted in the map: A) a spatial superspreader, 
which transmits over a large spatial area, potentially connecting sub-populations and may be important from an 
epidemiological perspective; and B) a superspreader with a limited dispersal range that infects a large number of progeny 
but remains within a small spatial radius. Inset map shows the location of the Serengeti District within Tanzania. 



prevention and control," August 2010, http://www. 
rabiesblueprint.com/, and Lembo et al. 2011). The 
rabies blueprint presents a major step in the global 
fight against rabies, providing guidelines and econ- 
omically feasible strategies to aid policy makers and 
local communities seeking to embark on rabies 
intervention and control measures. Landscape gen- 
etics research constitutes an important resource 
within this multidisciplinary approach to help advise 
and sustain successful control initiatives. 

Vaccination is widely deemed the most effective 
means of rabies control with demonstrated successes 
in reducing incidence and eliminating disease even in 
areas with limited resources (Cleaveland et al. 2003; 
Schneider et al. 2007; Lembo et al. 2010). Molecular 
genetics has demonstrated that domestic dogs are 
critical reservoirs for canine rabies, even in parts of 
Africa with abundant wildlife populations, and 
therefore indicate that controlling dog rabies 
through vaccination should eliminate infection from 
all other species (Lembo et al. 2007, 2008). By 
changing the susceptibility of populations, mass 
vaccinations change the landscape in which rabies 



circulates, and at sufficiently high levels of vacci- 
nation coverage, transmission can be interrupted. 
There are two mechanisms by which vaccination 
alters the landscape to control rabies. The first is that 
vaccination itself creates a barrier of susceptible 
individuals that block the dispersal of the pathogen 
across the landscape, and the second is via the 
reinforcement of existing barriers with vaccination. 
Using landscape genetics to explore the potential 
impacts of these aspects of landscape control is a 
logical next step. 

Recurrent rabies epidemics occur across large areas 
where vaccination programmes are patchy and 
unsustained (Hampson et al. 2007), suggesting that 
a proactive long-term vaccination programme 
coordinated across political boundaries is required 
for success. Indeed, the effectiveness of such a 
programme has been proven by the intensive 
vaccinations coordinated by the Pan American 
Health Organization (PAHO) in Latin America 
over the past few decades: dog rabies has been 
eliminated in a large portion of the southern 
continent, and reported cases from other countries 
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are highly localised due to restrictions on trans- 
mission pathways from vaccination barriers 
(Schneider et al. 2007). Crucially, cases of human 
rabies have dropped in these areas, reinforcing the 
importance of effective dog rabies control strategies. 
Strategically placing a vaccine barrier could maintain 
freedom from rabies resulting from successful inter- 
ventions in otherwise landlocked areas. 

Environmental or anthropogenic barriers to natu- 
ral transmission of rabies offer the opportunity to 
strengthen and smarten vaccination initiatives. Vacci- 
nation can be viewed as a form of barrier that impedes 
rabies spread, and therefore many of the techniques 
used to draw insights on the permeability of barriers 
could equally be applied to vaccination programmes. 
Oral Rabies Vaccination (ORV) campaigns for wild- 
life, including bait distribution in proximity to a 
pre-existing natural barrier (the Appalachian Moun- 
tains), and utilising existing environmental features 
to reinforce control campaigns (Wandeler et al. 1988) 
have successfully contained wildlife rabies. Past 
experience has shown vulnerability to breaches 
associated with the differential permeability of 
barriers and emphasizes the need for ongoing and 
targeted surveillance to enable early detection and 
swift responses to incursions (Russell et al. 2005). Yet 
despite this, few studies have quantified the utility 
of barriers within a landscape, and there are no 
guidelines available for the design and implemen- 
tation of effective cordons sanitaires for dog rabies. 

In addition to barrier studies, modern application 
of spatial data is allowing us to make increasingly 
accurate measurements of epidemiological par- 
ameters that may affect disease dynamics. For 
example, Bharti et al. (2011) demonstrate the use 
of remote sensing to test for human predictors of 
disease. They used anthropogenic light from satellite 
imagery as a measure of seasonal fluctuations of 
human populations. The observed fluctuation in 
light intensity (as a measure of population density) 
correlated to measles distribution and spread in cities 
in Niger, and provided an accurate, near real-time 
representation of short-term population fluctuations 
that may drive pathogen transmission. This approach 
demonstrates the use of relatively simple proxies for 
quantifying migration patterns in poorly resourced 
regions, often the same regions carrying the highest 
disease burden from dog rabies. 

Large-scale interventions are expensive, and adap- 
tive management is often required alongside inter- 
vention. Refinements in resource management will 
ultimately rely on a combination of knowledge 
from genetics, landscape and host ecology as part 
of a reactive programme. Specific landscape elements 
affecting dog rabies spread are generally not known 
a priori, so exploration of phylogeographic patterns 
may identify genetically distinct viral or host 
populations in areas, which could be targeted for 
vaccination. However, experience from wildlife 



rabies suggests that caution is warranted when taking 
such an approach. Firstly, apparent boundaries 
between areas dominated by different genetic lineages 
may have emerged during the initial invasion process 
and thus do not necessarily reflect areas of low 
permeability for the virus (Real et al. 20056; Biek and 
Real, 2010). Secondly, while the stability of these 
phylogeographic domains indicates a lack of mixing 
between them, this simply suggests that immigrating 
viruses find it difficult to invade areas with an already 
established focus, but does not signify the absence of 
viral immigration per se. Such areas with putatively 
self-contained endemic foci could therefore experi- 
ence a high risk of rabies re-emergence following 
successful eradication, unless vaccination effort 
remains high. Whether the second consideration 
equally applies to dog rabies is currently not clear as 
pertinent empirical studies are lacking. As with many 
research problems in landscape genetics, these types 
of question may be very productively studied using 
simulation tools (Epperson et al. 2010). Simulations 
may indicate the most effective location for vaccine 
corridors/barriers e.g. vaccination on the far side of a 
natural barrier (Russell et al. 2006), whereas meta- 
population models may elucidate areas with high 
connectivity that could be the source of persistence in 
endemic areas. Quantifying the degree of connec- 
tivity between sub-populations is of critical impor- 
tance to understanding how rabies is maintained 
across landscapes, and is an integral part of designing 
effective control interventions and cordons sanitaires 
that limit and contain the virus. 

Ongoing surveillance is crucial to the long-term 
success of rabies control, and considerable value can 
be added to surveillance initiatives through the 
incorporation of landscape genetics, with the poten- 
tial to determine the source of incursions and reveal 
transmission pathways. Epidemiological surveillance 
may facilitate active case detection and identification 
of circulating strains, helping to identify areas missed 
by vaccination. Using a molecular genetics approach 
during the 2007 FMD outbreak in the UK, allowed 
swift and effective containment of the outbreak 
by directing interventions to the critical areas 
(Cottam et al. 2008). For countries with limited 
resources, surveillance in areas with ongoing control 
programmes and retrospective analysis may be useful 
for identifying remaining foci of infection to be 
targeted by vaccination, sources of incursions, or 
spillover events into or from wildlife. In addition, 
metapopulation models may be utilized to predict the 
direction of spread of rabies in naive populations 
(Beyer et al. 2011), directing the location of sentinel 
points and control measures to 'hotspots' of infection 
(Haydon et al. 2006) and putative transmission 
networks based on sampled cases should indicate 
how well current levels of surveillance are capturing 
rabies incidence based on inferred missing links 
between sampled cases. 
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Rabies containment and control will require 
ongoing surveillance and sustained control efforts. 
In time, when rabies incidence is reduced to low 
levels (the period when control measures often lapse), 
genetics can provide a means of directing resources to 
where they are most needed, maintaining a cost- 
effective approach. Given the small genome of RABV 
(12Kb), and advances in NGS sequencing, future 
exploration of phylogeographic patterns utilizing 
whole genome sequencing is a realistic prospect. 
This promises to provide appropriately fine genetic 
resolution for samples collected on small spatio- 
temporal scales that might otherwise be uninforma- 
tive. Application of novel techniques in landscape 
genetics to other RNA viruses such as FMD 
and Influenza highlights interesting possibilities for 
uncovering disease dynamics and aiding control 
directives (Cottam et al. 2008, Bedford et al. 2010). 
The technology and analytical power are available 
but have yet to be fully exploited for dog rabies, and 
thus offer exciting prospects on what can be achieved 
with their implementation in the future. 

CONCLUSIONS 

There is evidently a wide scope for the use of 
landscape genetics to explore and understand the 
dynamics of pathogen spread and persistence. Dogs 
are the principal reservoir of rabies, responsible for 
the majority of human rabies cases, yet we know little 
about the dynamics of the pathogen in this host. In 
order to reduce the many thousands of human rabies 
deaths that occur annually due to contact with rabid 
dogs, it is crucial that we uncover the mechanisms 
governing viral dispersal across the landscape so as to 
direct successful control interventions. The foun- 
dation of knowledge from studies in wildlife popu- 
lations provides a starting point but aspects of host 
behaviour, the inherent influence of humans, and the 
long-term endemic nature of dog rabies foci requires 
a modified approach be taken for dog rabies. There 
lies great potential for advancing the effectiveness of 
control campaigns in areas burdened with disease; 
specifically landscape genetics has most to contribute 
to improving surveillance and modifying control 
strategies based on information gleaned from sur- 
veillance. This includes the design and placement of 
cordons sanitaires, the prioritization of effort towards 
persistent foci or targeting sources of outbreaks and 
conduits of transmission. Given the availability of 
powerful new genetic and spatial techniques, efforts 
now need to push towards real-world application of 
landscape genetics to rabies control and elimination. 
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